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We consider here the spin independent neutralino-proton cross section for a variety of SUGRA and 
D-brane models with R-parity invariance. The minimum cross section generally is ~ 1 X 10" -(9-l0)pb 
(and hence accessible to future detectors) except for special regions of parameter space where it may 
drop to ~ 10~ 12 pb. In the latter case the gluino and squarks will be heavy (~1 TeV). 



Dark matter detectors have now achieved 
a sensitivity that they have begun to probe 
interesting parts of SUSY parameter space. 
It is thus of interest to see what sensitivity 
will be needed to explore the entire space. To 
examine this we consider here models based 
on gravity mediated supergravity (SUGRA), 
where the LSP is generally the lightest neu- 
tralino (x?)' Neutralinos in the halo of the 
Milky Way might be directly detected by 
their scattering by terrestrial nuclear targets. 
Such scattering has a spin independent and a 
spin dependent part. For heavy nuclear tar- 
gets the former dominates, and it is possible 
to extract then (to a good approximation) the 
neutralino -proton cross section, <r^o_ p . Cur- 
rent experiments (DAMA, CDMS, UKMDC) 
have sensitivity to halo Xi for 

a^,_ p ~ 1 x lCn 6 pb (1) 

and future detectors (GENIUS, Cryoarray) 
plan to achieve a sensitivity of 

<Jyt{-v ~ ( 1(r9 ~ 10~ 10 )pb (2) 

We consider here two questions: (1) what 
part of the parameter space is being tested by 
current detectors, and (2) what is the small- 
est value of cr^o_ p the theory is predicting 
(i.e. how sensitive must detectors be to cover 
the full SUSY parameter space). The an- 
swer to these questions depends in part on 
the SUGRA model one is considering and 
also on what range of theoretical and in- 
put parameter one assumes. In the follow- 
ing, we examine three models that have been 



considered in the literature based on grand 
unification of the gauge coupling constants 
at Mq = 2 x 10 16 GeV: (1) Minimal su- 
per gravity GUT (mSUGRA) with universal 
soft breaking at Mq EJ; (2) Nonuniversal soft 
breaking models for Higgs and third genera- 
tion scalar masses at Mq, and D-brane mod- 
els (based on type IIB orientifoldsa) which 
allow for nonuniversal gaugino masses and 
nonuniversal scalar masses at Mq a. 

While each of the above models contain 
a number of unknown parameters, theories 
of this type can still make relevant predic- 
tions for two reasons: (i) they allow for ra- 
diative breaking of SU(2) x U(l) at the elec- 
troweak scale (giving a natural explanation 
of the Higgs mechanism) , and (ii) along with 
calculating a%o_ p , the theory can calculate 
the relic density of 

Xi, i-e 

where p^o is the relic mass density of Xi an d 
p c = 3Hq /8ttGn (Ho is the Hubble constant 
and Gat is the Newton constant). Both of 
these greatly restrict the parameter space. In 
general one has Q^oh 2 ~ (J Q Xf dx{a ann v))^ 1 
(where cr ann is the neutralino annihilation 
cross section in the early universe, v is the rel- 
ative velocity, Xf = kTf /m^o , Tf is the freeze 
out temperature, (...) means thermal average 
and h = Ho/100 km s^Mpc" 1 ). The fact 
that these conditions can be satisfied for rea- 
sonable parts of the SUSY parameter space 
represents a significant success of the SUGRA 
models. 

In the following we will assume Ho = 
(70 ± 10)km s _1 Mpc _1 and matter (m) and 
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baryonic (b) relic densities of fl m = 0.3 ± 0.1 
and n b = 0.05. Thus n^ah 2 = 0.12 ± 0.05. 
The calculations given below allow for a 2a 
spread, i.e. we take 0.02 < il^oh 2 < 0.25. It 
is clear that when the MAP and Planck sat- 
telites determine the cosmological parame- 
ters accurately, the SUGRA dark matter pre- 
dictions will be greatly sharpened. 



1 Calculational Details 

In order to get reasonably accurate results, it 
is necessary to include a number of theoret- 
ical corrections in the analysis. We list here 
the main ones used in the calculations be- 
low: (i) In relating the theory at Mq to phe- 
nomena at the electroweak scale, the two loop 
gauge and one loop Yukawa renormalization 
group equations (RGE) are used, iterating to 
get a consistent SUSY spectrum, (ii) QCD 
RGE corrections are further included below 
the SUSY breaking scale for contributions in- 
volving light quarks, (iii) A careful analysis 
of the light Higgs mass mh is necessary (in- 
cluding two loop and pole mass corrections) 
as the current LEP limits impact sensitively 
on the relic density analysis for tan/3 < 5. 
(iv) L-R mixing terms are included in the 
sfermion (mass) 2 matrices since they produce 
important effects for large tan/3 in the third 
generation, (v) One loop corrections are in- 
cluded to rrif, and m T which are again impor- 
tant for large tan/3, (vi) The experimental 
bounds on the b — > 57 decay put significant 
constraints on the SUSY parameter space 
and theoretical calculations here include the 
leading order (LO) and approximate NLO 
corrections. We have not in the following im- 
posed b — r (or t — b — r) Yukawa unification or 
proton decay constraints as these depend sen- 
sitively on unknown post- GUT physics. For 
example, such constraints do not naturally 
occur in the string models where SU (5) (or 
5O(10)) gauge symmetry is broken by Wilson 
lines at Mq (even though grand unification of 
the gauge coupling constants at Mq for such 



string models is still required). 

All the above corrections are under theo- 
retical control except for the b sj analysis 
where a full NLO calculations has not been 
done. (We expect that while the full anal- 
ysis might modify the regions of parameter 
space excluded by the b — » S7 experimen- 
tal constraint, the minimum and maximum 
values of cr% -p would probably not be sig- 
nificantly changed.) The analysis of a^o_ p , 
taking into account the above theoretical cor- 
rections has now been carried out by several 
roups obtain g results in general agreement 
These results are presented be- 
low. 

Accelerator bounds significantly limit the 
SUSY parameter space. Ir± the following we 
assume the LEP bounds @ m/, > 104(100) 
GeV for tan/3=3(5) and to + > 104 100 
GeV. (For tan/3 > 5, the bounds do 
not produce a significant constraint!!! ) For 
b — > S7 we assume aiL allowed range of 2a 
from the CLEO data t 2 ! The Tevatron gives 
a bound of m~ g > 270 GeV( for = m,g)H. 

Theory allows one to calculate the Xi~ 
quark cross section and we follow the analysis 
of E-3 to convert this to Xi ~ P scattering. For 
this one needs the ir — N a term, a^N and 
Co = CnN — (did + md){p\ss\p) and the quark 
mass ratio r = m s /(l/2)(m u + mj). We use 
Jiere a^N = 65 MeV, from recent analyses 
E§li3 based oil new n — N scattering data, 
a a = 30 MevQ and r= 24.4 ± 1.50. 



2 mSUGRA model 

We consider first the mSUGRA model where 
the most complete analysis has been done. 
mSUGRA depends on four parameters and 
one sign: mo (universal scalar mass at Mq), 
TO1/2 (universal gaugino mass at Mq), Aq 
(universal cubic soft breaking mass) , tan/3 = 
{H%)/{H\) (where (i/2,1) gives rise to (up, 
down) quark masses) and /V|/x|( where /1 is 
the Higgs mixing parameter in the superpo- 
tential, — fiHiH 2 ). One conventionally 



restricts the range of these parameters by 
"naturalness" conditions and in the follow- 
ing we assume m < 1 TeV, mi 12 < 600 GeV 
(corresponding to rrig < 1.5 TeV, m%a < 240 
GeV), \A /m \ < 5, and 2< tan/3 < 50. 
Large tan/3 is of interest since SO(10) mod- 
els imply tan/3 > 40 and also <r^o_ p increases 
with tan/3. &x a - P decreases with m^2 for 
large m^, and thus if one were to increase 
the bound on m 1/2 to 1 TeV (m s < 2.5 TeV), 
the cross section would drop by a factor of 2- 
3. 

The maximum o^o_ p arise then for large 
tan/3 and small mj^. This can be seen in 
Fig.l where (a yO ) max 

is plotted vs. rn^o 
for tan/3=20, 30, 40 and 50. Current detec- 
tors obeying Eq (1) are then sampling the 
parameter space for large tan /3, small m^o 
(and also small f^o/1 2 ) i.e 

tan/3 ~ 25; m^o - 90 GeV; fl^oh 2 ~ 0.1. 

(3) 

To discuss the minimum cross section, it 
is convenient to consider first m^o ~ 150 GeV 
( m i/2 < 350) where no coannihilation occurs. 
The minimum cross section occurs for small 
tan/3. One finds 

a^o_ p ~4x 10~ 9 pb; m^o ~ 140GeV (4) 

which would be accessible to detectors that 
are currently being planned (e.g. GENIUS). 

For larger m^o, i.e. mi/2 ~ 150 the phe- 
nomena of coannihilation can occur in the 
relic density analysis since the light stau, ?i, 
(and also e^, flu) can become degenerate 
with the Xi- The relic density constraint can 
then be satisfied in narrow corridor of mo of 
width Arno ~ 25 GeV, the value of mo in- 
creasing as mi/2 increases □. Since mo and 
mi/2 increase as one progresses up the corri- 
dor, cr^o „ will generally decrease. 

We consider first the case /i > 0LU. One 
finds in general that 0^0 _ p also decreases as 
A increases. Fig. 2 shows <J^o_ p in the do- 
main of large Aq and for two values of tan/3. 
One sees that the smaller tan/3 still gives the 
lower cross section, though the difference is 



mostly neutralized at larger mi/ 2 - (For large 
tan/3, m also becomes large to satisfy the 
relic density constraint i.e m = 700 GeV for 
tan/3=40, m 1 / 2 = 600 GeV.) We have in gen- 
eral for this regime 

a^o_ p ~ 1 x 10 _9 pb; form 1/2 < 600GeV, 
A« > , A a < 4m 1/2 . (5) 

This is still within the sensitivity range of 
proposed detectors. 

When fi is negative an "accidental" can- 
cellation can occur in part of the parameter 
space in the coannihilation region which can 
greatly reduce cr%°-p 0. This can be seen 
in Fig. 3, where starting with small tan/3 the 
cross section decreases, leading to a mini- 
mum at about tan/3=10, and then increases 
again for larger tan/3. At the minimum one 
has cTgo—p — 1 x 10~ 12 when tan/3=10 and 
TO i/2 = 600 GeV. More generally one has 

a%o_ p < 1 x 10~ 10 pb (6) 

for the parameter domain when 4 ~ tan/3 ~ 
20, m 1/2 ~ 450GeV(m§ ~ l.lTeV), /i < 0. 
In this domain, a^o_„ would not be accessi- 
ble to any of the currently planned detectors. 
However, mSUGRA also then predicts that 
this could happen only when the gluino and 
squarks have masses greater than 1 TeV (and 
for only a restricted region of tan/3) a result 
that could be verified at the LHC. 

3 Nonuniversal SUGRA Models 

In the discussion of SUGRA models with 
nonuniversal soft breaking, universality for 
the first two generations of squark and slep- 
ton masses at Mq is usually maintained to 
suppress flavor changing neutral currents. 
One allows, however, the Higgs and third gen- 
eration squark and slepton masses to become 
nonuniversal. We maintain gauge and gaug- 
ino mass unification at Mq- 

While these models contain a large num- 
ber of new parameters, their effects on cr^o_ p 
can be charcterized approximately by the 
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signs of the deviations from universality 
One choice can greatly increase cr^o_ p , by a 
factor of 10-100 compared to the universal 
case, and the reverse choice can reduce ct^o_,„ 
( though by a much lesser amount). Thus 
it is possible for detectors to probe regions 
of smaller tan/3 with nonuniversal breaking, 
and detectors obeying Eq. (1) can probe part 
of the parameter space for tan/3 as low as 
tan (3 ~ 4. 

The minimum cross section occurs (as 
in mSUGRA) at the lowest tan/3 and at the 
largest mi/ 2 i.e. in the coannihilation region. 
We limit ourselves here to the case where 
only the Higgs masses are nonuniversal. One 
finds then results similar to mSUGRA i.e. 
er^o_ p ~ 10~ 9 pb for fi > 0, < 600 

GeV. For a < 0, there can again be a cancel- 
lation of matrix elements reducing the cross 
section to 10~ 12 pb when m\i% — 600 GeV in 
a restricted part of the parameter space when 
tan/3 ~10. 

4 Summary 

We have examined here the neutralino-proton 
cross section for a number of SUGRA type 
models. In all the models considered, there 
are regions of parameter space with \i ~ P 
cross sections of the size that could be ob- 
served with current detectors. Thus with the 
sensitivity of Eq. (1), detectors would be 
sampling regions of the parameter space for 
mSUGRA where tan/3 ~ 25, m^o ~ 90GeV 

and Q^oh 2 ~ 0.1. Nonuniversal models can 
have larger cross sections and so detectors 
could sample down to tan/3 ~ 4, while for the 
D-brane models considered, detectors could 
sample down to tan/3 ^15. 

The minimum cross sections these mod- 
els predict are considerably below current 
sensitivity. Thus for mSUGRA one finds for 
a > that c-£0_ p ~ 1 x 10~ 9 pb for m 1( / 2 < 
600GeV, fi > 0, where m 1/2 = 600 GeV cor- 
responds to m g = 1.5 TeV, m^o = 240 TeV. 
This is still in the range that would be ac- 
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cessible to detectors being planned (such as 
GENIUS or Cryoarray). For u < 0, a cancel- 
lation can occur in certain regions of param- 
eter space allowing the cross sections to fall 
below this. Thus 

a^o_ p < 1 x 10~ 10 pb for 4 ~ tan/3 ~ 20, 
H < 0, m 1/2 ~ 450 GeV 

and reaching a minimum of a^o _ p = 1x10 
pb for tan/3=10, m 1/2 =600 GeV, [X < 0. This 
domain would appear not to be accessible to 
future planned detectors. Since m 1 / 2 =450 
GeV corresponds to rrig = 1.1 TeV, this re- 
gion of parameter space would imply a gluino 
squark spectrum at the LHC above 1 TeV. 

The above results holds for the mSUGRA 
model. While a full analysis of coannihila- 
tion has not been carried out for the nonuni- 
versal and D-brane models, results similar to 
the above hold for these over large regions 
of parameter space. Thus for nonuniversal 
Higgs masses and for the D-brane model one 



finds CT-£o_ p ~ 10 9 pb for a > 0, while a 
cancellelation allows a^a_„ to fall to 10~ 12 
pb for u < at tan/3 ~ 10 (with again a 
gluino/squark mass spectrum in the TeV do- 
main) . 
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Figure 1. (ct s o_ p ) 



max for mSUGRA obtained by 
varying Ao over the parameter space for tan/3 = 20, 
30, 40, and 50[9]. The relic density constraint has 
been imposed. 




Figure 2. (&^o_ p ) for mSUGRA in the coannihila- 
tion region for tan/3 = 40 (upper curve) and tan/3 = 3 
(lower curve), Aq = 4m 1 /2i /* > 0. 
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Figure 3. {a-o_ p ) for mSUGRA and fi < for 
(from top to bottom on right) tan/3=20, 5 and 10. 
Note that for tan/3 > 10, the curves terminate at the 
left due to the b — > S7 constraint. 



